On the existence of polynomial first integrals of 
quadratic homogeneous systems of ordinary differential 

equations 

Tsygvintsev Alexei 

Abstract 

We consider systems of ordinary differential equations with quadratic homogeneous right hand side. 
We give a new simple proof of a result already obtained in [8,10] which gives the necessary conditions 
for the existence of polynomial first integrals. The necessary conditions for the existence of a polynomial 
symmetry field are given. It is proved that an arbitrary homogeneous first integral of a given degree is a 
■ linear combination of a fixed set of polynomials. 
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1 Introduction 

In the present paper we study the system of ordinary differential equations with quadratic homogeneous 
right hand side 



X{ fi{p£\i • • • , 2-n)i fi ^ Q'ijk-Kj'Eki ^ijk G ^ 1, . . . , Tt. (1.1) 

Systems of such a form arise in many problems of classical mechanics: Euler-Poincare equations on Lie 
algebras, the Lotka-Volterra systems, etc. 

The main concern of this paper is to find the values of the parameters a^-fe for which equations (1.1) have 
first integrals. 

In the paper [2] the necessary conditions are found for the existence of polynomial first integrals of the 
system 

±i = Vi(xi, . . . ,x n ), i = l,...,n, (1.2) 

where Vi G C[xi, . . . ,x n ] are homogeneous polynomials of weighted degree s G N. In the case s = 2 we 
obtain equations (1.1). 

The method given in [2] is based on ideas of Darboux [1,6,7] who used a special type of particular solutions 
of the system (1.2) 

Xi{t) = di(/>(t), i = 1, . . . ,n, 

where <j>(t) satisfies the differential equation cf> = Atfi s , A is an arbitrary number and d = {d\, . . . , d n ) T ^ is 
a solution of the following algebraic system 



Vi(d) = Xdi, i = 1, . . . ,n. 



J-j | In this paper we generalize this method. 

It was shown in [8,10] that the weighted degree of a polynomial first integral of the system (1.1) is a 
certain integer linear combination of Kovalevskaya exponents (see [9]). In Section 2 we give a new simple 
proof of this result. In Section 3 a similar theorem for polynomial symmetry fields is proved. As an example, 
we consider the well known Halphen equations. Section 4 contains our main result. We present so called 
base functions and prove that every homogeneous polynomial first integral of a fixed degree is a certain 
linear combination of the corresponding base functions. In Section 5 we give an application of previous 
results to planar homogeneous quadratic systems where necessary and sufficient conditions for the existence 
of polynomial first integrals in terms of Kovalevskaya exponents are found. 

2 The existence of polynomial first integral. Necessary conditions 

Following the paper [2] , we consider solution C — • ■ ■ ? c n ) T ^ (0, . . . , 0) T of algebraic equations 

fi(ci,...,c„) +Ci = 0, i=l,...,n. (2.1) 
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Kij = — -(C) +Sij, i,j = l,...,n. 



Define the Kovalevskaya matrix K [3] 

dxj' 

where Sij is the Kronecker symbol. Let us assume that K can be transformed to diagonal form 

K = diag(pi, . . .,p n ). 
The eigenvalues pi , . . . , p n are called Kovalevskaya exponents. 

Lemma 1. ([3]) Vector C is an eigenvector of the matrix K with eigenvalue p\ = —1. 
Consider the following linear differential operators 

i—1 1 i->3 = 1 1 

i=i Ud - 1 i=i UXt 



(2.2) 



which satisfy relations 

[D-,D+] = D -U, [D ,D-]=D-, (2.3) 

where [A, B] = AB - BA. [8,10] 

Theorem 2. Suppose that the system (1.1) possesses a homogeneous polynomial first integral Fm of degree 
M, and pi = —l,p2,...,p n are Kovalevskaya exponents. Then there exists a set of non-negative integers 
tc2, ■ ■ ■ , k n such that 

k 2P2 + ■■■ + k nPn = M, k 2 + --- + k n <M. (2.4) 
Proof By definition of a first integral D + Fm = 0. Considering identities 

DL(D+F m ) = 0, for ZeN, (2.5) 
we obtain the following set of polynomials 

Fm,Fm-i, • ■ • , F p+ i,F p , 

defined by the recursive relations 

D-F i+1 = {M-i)F u i = p,...,M-l, 

where the number 1 < p < M is determined by the condition D-F p = 0. Using (2.3), (2.5) we deduce the 
following chain of equations 

DqF m = MF M - D+Fm-i, 
DqFm-i = MFm-i - D+Fm-2, 

D F p+1 =MF p+1 -D+F p , 

D F p - MF p . (2.6) 

Let Ji , ... , J n - linearly independent eigenvectors of the Kovalevskaya matrix K corresponding to the 
eigenvalues p\ = — 1, p 2 , ■ ■ ■ , p n - According to Lemma 1 we can always put J\ = C. 
Consider the linear change of variables 

n 

Xi = ^L ij y j , i = l,...,n, (2.7) 

3=1 

where L = (Lij) is a nonsingular matrix defined by 

L = (C, J 2 , ■ ■ ■ , Jn), 
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then obviously 

L~ X KL = diag(-l, p 2 , . ■ . , p„). 

With help of (2.7) and Lemma 1 one finds the following expressions for the operators D , D- in the new 
variables 

" d d 

D o = y2piVt^-, D - = -^—> 



and the equation (2.6) becomes 



(*m-n- + --- + PnV n -n-)F p = MF p . (2.8) 
oy 2 oy n J 

We can write the polynomial F p as follows 

F p = ]T A k2 ... kn y k ^ ■■■y k n ^ \k\ = k 2 + ■ ■ ■ + k n , fc,eZ + . (2.9) 

\h\= P 

Substituting (2.9) into (2.8), one obtains the following linear system 

(k 2 p 2 + ■■■ + k nPn )A k2 ... kn = MA k2 ... kn , for \k\= p. (2.10) 

Taking into account that F p is not zero identically, we conclude that there exists at least one nonzero set 
k 2 , . . . , k n , \k\ < M such that 

k 2 p 2 + ■■■ + k nPn = M. (2.11) 

This relation implies (2.4). Q.E.D. 

Remark. Theorem 2 does not impose any restrictions on grad(F^) calculated at the point C. Thus, it 
generalizes the theorem of Yoshida ([3], p. 572), who used essentially the condition grad(F^) ^ 0. 



Corollary 3. The Halphen equations 



x\ = x 3 x 2 - xxx 3 - xix 2 , 

x 2 = xix 3 - x 2 x\ - x 2 x 3 , (2-12) 
x 3 = x 2 x\ - x 3 x 2 - x 3 x\, 



admit no polynomial first integrals. 



Indeed, the system (2.12) has Kovalevskaya exponents p\ = p 2 = p 3 = — 1. It is easy to verify that 
conditions (2.4) arc not fulfilled for any positive integer M. Moreover, as proved in [2], the system (2.12) 
has no rational first integrals. 



3 Existence of polynomial symmetry fields. Necessary conditions 

The first integrals are the simplest tensor invariants of the system (1.1). In [4] Kozlov considered tensor 
invariants of weight-homogeneous differential equations which include the system (1.1). In particular, he 
found necessary conditions for the existence of symmetry fields. Below we propose a generalization of his 
result. 

Recall that the linear operator W — w i( x i, ■ ■ ■ , x n)j¥—, is called the symmetry field of (1.1), if 

i=i OXi 

[W, D + ] = 0, where D + is defined by (2.2). If w\, . . . ,w n are homogeneous functions of degree M + 1 then 
the degree of W is M ([4]). 

Theorem 4. Suppose that the system (1.1) possesses a polynomial symmetry field of degree M and 

pi = —l,p 2 ,...,p n are Kovalevskaya exponents. Then there exist non-negative integers k 2 ,...,k n , \k\ < 

M + 1 such that at least one of the following equalities holds 

k 2 p 2 -\ \-k n p n = M + pi, i=l,...,n. (3.1) 
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Proof. Let Wm be a polynomial symmetry field of degree M. Substituting in the proof of Theorem 
2 operators D + , D , D- with their commutators [-D+,], [Do,], [D-,] respectively we repeat the same 
arguments. 

Corollary 5. The Halphen equations (2.12) admit no polynomial symmetry fields. 

Using (3.1) we obtain that M may be equal to — 1, 0, 1 only. It is easy to check that (2.12) does not 
have symmetry fields of such degrees. 



4 Base functions 

After the change of variables (2.7) the system (1.1) takes the form 

yi = -y\ + <pi(y2,---,y n ) , 41 ^ 
Vi = (pi - l)yiyi + <Pi(y 2 ,-- -,y n ), i = 2,...,n, 

where ipi are quadratic homogeneous polynomials in the variables y 2 , • • • , y n - 
According to (2.2) define operators D+, Do, D-. 

A homogeneous polynomial Pm(vi, ■ ■ ■ , Vn) of degree M satisfying the condition 

D_(D+P M )=0, (4.2) 

is called base function of the system (4.1). In others words, the function D + Pm does not depend on y\. It 
is clear that base functions of degree M form a linear space Lm over field C. 

Lemma 6. If the system (4-1) has a homogeneous polynomial first integral Fm of degree M , then Fm G Lm- 
Indeed, by definition, we have D + Fm — 0, hence, in view of (4.2), Fm G Lm- 

Let J(M) = {z G Z" _1 I z 2 p2 + • • • + z n p n = M, \z\ < M} be the set of integer-valued vectors z = 
(z 2 , . . . , z n ) T for which the condition (2.11) is fulfilled. Put m = \ J(M)\ and suppose J(M) ^ 0. 

Theorem 7. The dimension d of Lm satisfies the condition 1 < d < m. 

Proof. Let us assume the set J(M) contains vectors z^\ . . . , z^ which are ordered by the norm 

\z\ = z 2 H V z n 

\ Z W\ < ... < \ z (m )\. 

Define the vector p — (p 2 , ■ ■ ■ , p n ) T and put (p, z) = p 2 z 2 + • • • + p n z n , \ = rij, i — 1, . . . , m. 

Following the proof of Theorem 2, for each i = 1, . . . , m consider the system of linear partial differential 
equations 

DoPi, ni =MP i!ni , 

DoP-L ni + l = MPi, ni + i — D + P iyUi , 



DoPi,M-l — MPi t M-l — D + P iy M- 2 , 
DoPi,M — MP iy M — D + Pi t M-l, 



(4.3) 



D_P M+1 = (M - l)P itl , l = m,...,M-l, (4.4) 

which defines polynomials Pj i7li , . . . , P^m recurrently. 

It follows from (z [% > , p) = M that the first equation in (4.3) has the particular solution Pi, ni = y 2 • • • yn ■ 
Equations (4.3), (4.4) define certain base function P, ; m- Indeed, according to (4.4), we have 

P t M-i - D_P iM . (4.5) 
Substituting (4.5) into the last equation in (4.3), and using relations (2.3) we find 

D Pi, M = MP hM - D+D-P hM = MP lM - (D_D+ - D + U)P lM - 
Hence D_(D + P iiM ) = 0. 
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Now consider the problem on the existence of a solution of (4.3), (4.4) in form of homogeneous polynomials 
Pi iTli , . . . , Pi,M- Fix certain i = 1, . . . , m and put Oj = M — n$. Using the relations (4.4) we can write 

p ^\ (p—j\ p—j j -i (4-6) 

*i,ni+p — \ a -j)Vl l-i.nt+j, P~ U---? a i? 

j=0 1 

where Ii t k{V2, ■ ■ ■ , Vn) are certain homogeneous polynomials of degrees k = rij, . . . , M. Notice that I^k does 
not depend on y\. 

Differential operators D +1 Dq can be represented in the form 



D + = {-y\ + Vl )JL + yi (A -U) + A+, 



(4.7) 



m 

where 

^+ = E^> A ° = £^^ ^ = 5>aw ( 4 - 

fe=2 WK fe=2 WK fc=2 yK 

Using (4.3), (4.6), (4.7) one deduces the following equations for determination of I 



A Ii, nt = MI i<ni , 

AoIi,m+2 = MIi yTli+ 2 — a.i<p\Ii,m — A.+i%,m+\i 
AoIi,m+3 = MIi,n t +3 - (a t - l)<piJi, ni +i - A + I itni+2 , 



(4.9) 



AqI^m — MI iy M — 2<PlIi,M-2 — A + Ii t M-l- 

We can write each equation of (4.9) as follows 

AoXi=MXi + Yi, (4.10) 
where Xi, Y; are homogeneous polynomials of weighted degree I = m, . . . , M. Let us assume 

Xi = ^2 c i2 ... in yi 2 ■■■y%, Yt = d i2-i n y l 2 2 1*1 = *2 + •■• + *«, (4.11) 

|i|=/ |i|=; 

where Cj 2 ...j n , e?i 2 ...i„ are constant parameters. Then substituting (4.11) into (4.10), we obtain the following 
linear system with respect to Cj 2 ...j n 

(«2P2 H - M)c i2 ...j„ = di 2 ... in , (4.12) 

for i 2 , . . . ,i n = 0, 1, . . . , |i| = I. 

Suppose there exists a set k 2 , ■ ■ ■ , k n for which the following conditions are fulfilled 

fe...,y T £J(M), d k2 ... kn ^0, \k\=l, (4.13) 

then the solution Ii. ni , • • • , h,M does not exist. In this case we put P^m = 0. 
If the conditions (4.13) are not satisfied, we obtain the base function 

at 

I'..m X ' / < " - * '-I 1' 

3=0 

It is easy to show that polynomials {Pi,M}lZm are linearly independent over the field C. 

Taking into account that ni < • • • < n m and using (4.13), we see that in case i = m we always can 
determine the base function P^m- Therefore, under the assumption J(M) ^ 0, the space Lm always 
contains a nonzero function. Q.E.D. 

Corollary 8. If at least one resonance condition of the form 

(z,p)=M, \z\<M, zeZl~\ 
is fulfilled, then there exists a base function of degree M. 
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5 Polynomial first integrals in the case of quadratic plane vector 
field. 

The first classification of integral curves of two-dimensional quadratic homogeneous systems can be found in 
the paper of Lyagina [5] and later was completed by numerous authors. 

In this section we apply the previous results to this problem to illustrate the method of basis functions. 

Consider the system 

x\ = a\x\ + b\X\X2 + d?,x\ 

±2 = 0,2X\ + biX\X2 + dixl, 

where ai,bi,di are constant parameters. 

Let = (ci~\c^) T , = {c"i\c 2 2 " > ) T be any two linearly independent solutions of the algebraic 
system (2.1). The exceptional cases when the system (2.1) has only one or admit no solutions are excluded 
for the discussion below. 

Assume that Kovalevskaya exponents corresponding to , are 

- (-l,Pi) T , 3*2 = (-l,p 2 ) T (5.2) 

Lemma 9. The system (5.1) has a homogeneous polynomial first integral of degree M if and only if there 
exists integer k = 1, . . . , M — 1 such that pi = M jk and p2 is one of the following numbers 

MM M 

M. 



M-fc'M-fc-1''"' 2 
Proof. Consider the following change of coordinates 



J \ 4 1} c 2 ] I y l>2 



(5.3) 



2 , 1W „. (5-4) 



which exists because of linear independence of vectors . In coordinates (pi,P2) the system (5.1) 

takes a more simple form 

Pi = ~p\ + (P2 - l)PlP2, 

P2 = -P\ + (Pl - l)PlP2- 

It is easy to show that under the change (5.3), the vectors turn into = (1,0) T , = (0, 1) T 

respectively. Obviously, the system (5.4) has the same Kovalevskaya exponents (5.2). The matrix K, 
calculated for is 

' -1 02 -1 
P i 



K 



Under the assumption pi ^ —1, we can reduce K to a, diagonal form using the following change of 
coordinates 

' Pl ) = L ( yi 

Pi V 2/2 



with the constant matrix L 



The case pi = — 1 will be considered below. 
Finally, equations (5.4) take the form (4.1) 



1 P2-1 
pi-l 



V2 = (Pl ~ 1) J/1 2/2 + V>2, 



where 

fi=ay2, W2 = byl, 

a = (P2 - l)(pi + P2), b= (pi - l)(p 2 - 1) - pi - 1. 
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For the operators (4.8) we get 

9 9 T~T d 

A+ = (fi2-K—, A = p 1 y 2 — , U = y 2 — . 
oy 2 oy 2 oy 2 

Let Fm be a polynomial first integral of (5.5) of degree M. 

According to Theorem 2, there exists an integer k = 1, . . . , M — 1 such that 

k Pl = M. (5.6) 

We exclude the case k = M(p\ = 1), since if pi = ±1, then the system (5.1) has no polynomial first integrals. 
This can be shown directly using equations (5.4), (5.5). 

Next, calculate the base function Pm corresponding to the resonance condition (5.6). Consider the 
equations (4.9). It is obvious, that polynomials . . . , J^m can be represented in the following form 

h,k+i = aiy 2 +i , i = 0,...,M-k, (5.7) 

where ceo, ■ ■ ■ , &M-k are constant parameters. 
Substituting (5.7) into (4.9) we obtain 

a a = 1, 

bk 

ai = pi (k + 1) — M ' (5.8) 

_ a(M -k-i + 2)a t _i + bjk + i - l)q t _ 2 , 
a.- p-^k+Iy^M ' * k - 

According to (4.14), we get the following expression for the base function Pm 



M-k 

°M = 



y: a jy r k - j y k 2 +j - (5.9) 

3=0 



By definition of the base function it is clear that 

D+P M = 6y? +1 , (5.10) 

where 

5 = actM-i + bMoeM- (5-11) 

Thus, the linear space Lm contains only one polynomial Pm- Hence, taking into account Lemma 6, Fm = 
const • Pm- 

Using (5.10), we conclude that Pm is a first integral if and only if S = 0. In view of (5.6), (5.8), (5.11) 
and the above condition, we arrive at Lemma 9. 

Theorem 10. The system (5.1) possesses a homogeneous polynomial first integral of degree M if and only 
if the following conditions are fulfilled 

a) pi, i — 1,2 - are positive rational numbers, 

Vp^ + pz 1 < i, 

c)^N. 
Pi 

This is an obvious consequence of Lemma 9. 
As an example consider the following system 



X ^ — X 9 X 2 f 

x 2 = — 3xf — 8xix 2 + 3x2- 



(5.12) 
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The vectors c^ x \ have the form 

c« = (1/8, -1/8) T , = (1/8,1/8) T . 

Calculating the corresponding Kovalevskaya exponents (5.2) one obtains 

$ft 1 = (-l,3) T , 5R 2 = (-1,3/2) T . 

We have p\ = 3, pi — 3/2, p^ 1 + p^ 1 = 1. So, the conditions a), b) of Theorem 10 arc fulfilled. By the 
condition c) one gets M — 31, I 6 N. Thus, the equations (5.12) possess a cubic first integral F 3 . Using 
formulas (5.8), (5.9), we obtain 

F3 — x i+ x\x2 — X\x\ — x\. 
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